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Abstract. - The Hamiltonian of a polariton model for an inhomogeneous linear absorp- 
tive dielectric is diagonalized by means of Fano's diagonalization method. The creation and 
annihilation operators for the independent normal modes are explicitly found as linear com- 
binations of the canonical operators. The coefficients in these combinations depend on the 
tensorial Green function that governs the propagation of electromagnetic waves through the 
dielectric. The time-dependent electromagnetic fields in the Heisenberg picture are given in 
terms of the diagonalizing operators. These results justify the phenomenological quantization 
of the electromagnetic field in an absorptive dielectric. 



Introduction. - To describe quantum optical phenomena in dielectrics it is essential 
to have available a quantization procedure for the electromagnetic field in ponderable mat- 
ter. Preferably, such a quantization should be based on the standard canonical quantization 
method of quantum field theory. To apply that method to the fields in linear dielectrics, 
one should start from a Hamiltonian description, in which both the field and the dielectric 
are given in terms of canonical variables that arc coupled bilincarly. The normal modes in 
such a system are the well-known polaritons [1]. If the dielectric is absorptive and disper- 
sive, damping of the polariton modes may be taken into account by coupling the dielectric 
degrees of freedom to a suitable bath of harmonic oscillators with a continuous range of fre- 
quencies. Working along these lines, Huttner and Barnett [2] were the first to formulate a 
damped-polariton model for an absorptive dielectric in an electromagnetic field and to study 
its properties. By using a diagonalization method due to Fano [3] they were able to find the 
full time dependence of the electromagnetic field operators for their model. 

The method of diagonalization employed in [2] is based on a separation of longitudinal 
and transverse degrees of freedom, and on a Fourier decomposition of the canonical variables. 
Both of these means are only expedient for homogeneous systems with translation invariance. 
In fact, when the damped-polariton model is taken to be inhomogeneous, with an arbitrary 
spatial dependence of the material properties, the longitudinal and transverse degrees of free- 
dom get coupled, while the Fourier components of the canonical variables for different wave 
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vectors start interacting as well. Hence, the diagonalization procedure in [2] runs into diffi- 
culties for the inhomogeneous case. Since inhomogencitics are unavoidable in any quantum 
optical experiment involving dielectrics, this is a serious drawback. It is the purpose of this 
Letter to show that diagonalization of the inhomogeneous version of the damped-polariton 
model is possible, and that explicit expressions for the diagonalizing operators in terms of 
canonical variables can be found. 

Model. - In the damped-polariton model the polarization density is coupled to a bath 
of harmonic oscillators with a continuous range of eigenfrequencies. The electromagnetic 
field interacts with the dielectric according to the standard minimal- coupling scheme. The 
Hamiltonian of the model is [2] 

/•OO 2 1 rOQ 

Jo P 2p pJq 

+ / drdr' ^- — ^- ^ . (1) 

J 87reo|r-r'| ^ ' 

The transverse part of the electromagnetic field is determined by the vector potential A(r), for 
which the Coulomb gauge is adopted. Its conjugate canonical momentum is n(r). The linear 
dielectric, with a space-dependent density p{r), is described by the harmonic displacement 
variable X{r) and its canonical momentum P{r). The associated (renormalizcd) eigenfre- 
quency LOo{r) is generally space-dependent as well. The electromagnetic field is coupled to 
the dielectric variable X in the usual way. In terms of the polarization density —aX, with 
a space-dependent coupling parameter a(r) > 0, the minimal coupling scheme leads to an 
electrostatic contribution and to a bilinear interaction term with A ■ P. Finally, damping is 
introduced in the system by a continuum bath of harmonic oscillators with canonical variables 
Y^{r), Quiir) and with eigenfrequencies uj. These bath oscillators are coupled to X{r) with 
a strength Vu{r) > 0. We used the notation X' = X{r'), and likewise a' and V. 
The canonical variables obey the standard commutation relations 

[U{r),A{r')] = -ihST{r-r'), [P{r),X{r')] = -ih\S{r - r') , 
[QAr),Y^'{r')] = -ih6{co-w')\5{r-r'), (2) 

while all other commutators of the canonical variables vanish. Here I is the three-dimensional 
unit tensor, while 6T{r) = I S{r) + W(47rr)~-'^ is the transverse delta bmction. 

The electric field operator E is the sum of a transverse part depending on 11 and a 
longitudinal part that is proportional to the polarization density: 

E{r) = -—U{r) + — [aX{r)]L. (3) 

So So 

The longitudinal part of a vector (or a tensor) is obtained by a convolution with the longitu- 
dinal delta function <5l('") = — VV(47rr)~^. The displacement field 

D(r) = -n(r) - [aX{r)]T (4) 

is purely transverse. 

In the following we will show how the Hamiltonian (1), with canonical operators satisfying 
the commutation relations (2), can be brought in diagonal form. 
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Fano diagonalization. - The Hamiltonian is quadratic in the canonical variables. Hence, 
it should be possible to find a diagonal representation of the form 



dr J dojhojC\r,uj)-C{r,uj), 



(5) 



where we omit a zcro-point-energy term. The operators C(r,w) are annihilation operators, 
which (together with the associated creation operators) satisfy the commutation relations: 

[C{r, Lo), Ct(r', w')] = ^('^ - I ^(^ - r') , [C{r, to), C{r', to')] = . (6) 

Each canonical operator can be written as a linear combination of the annihilation and 
creation operators. For instance, one has 

A{r) = J dwf^(r,r',w)-C(r',w)+h.c., 

QUr) = J dr' J^^ doj'fQ{r,r',oj,uj')-C{r',u')+h.c., (7) 

with tensorial coefficients f^ and fg. The coefficients fn, fx, fp and fy are defined analogously. 
Both f^ and fn are transverse in r. The electric field follows from (3) as 



/nOO 
dr' J dwfE(r,r',w) •C(r', w)+h.c., 



(8) 



with the coefficient 



fE{r,r',uj) = fn{r,r',oj) + — [afx{r,r' ,oj)]^ . (9) 

From eq. (6) it follows that the coefficients are equal to commutators. For instance, one 
has: 

f^(r,r',w) = [A(r),Ct(r',a;)] , fQ{r,r' = [QUr),C\r' ,cj')] . (10) 
Inversely, each C{r, ui) is a linear combination of the canonical operators: 

C{r,u,) = Jdr'S^A{r')-f^{r',r,i,)-U{r')-fX{r',r,u;) 
+X{r')-rp{r',r,Lu) - P{r') ■ r^ir' ,r,u;) 

+ jj^dw' [Y^,{r')-rQ{r',r,u',u)-Q^,{r')-f;.{r',r,iv',u)]Y (11) 

Fano's method to diagonalize the Hamiltonian amounts to finding the tensorial coefficients in 
these expressions by solving a set of equations that follows from the commutator of C with 
the Hamiltonian. In fact, eqs. (5) and (6) imply: 

[C{r,aj),H]=h(jC{r,uj). (12) 

Upon inserting eqs. (1) and (11), employing (2) and comparing the coefficients of the canonical 
operators, we arrive at the following set of linear relations: 

iwfA{r,r',u!) = — -fn{r,r',w) , (13) 
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iufx{r,r',u}) 
iujfp{r,r',Lo) 

io^fy (r, r', Lu' , ui) 
iu;fQ{r,r',u;',uj) 



AfA{r,r',oj) + 

Mo 



■fA{r,r',u;) 



■fp{r,r',Lo) 



--fA{r,r',ui) - -fp(r,r',a;), 

P P 

P^l^x{r,r' ,uj) A [afx(r,r', 

1 

P Jo 

— Wc^' fx(r,r',w) fQ(r,r', . 

pa;'^fy(r,r',a;',w) , 



(14) 
(15) 

(16) 

(17) 

(18) 



with a, p, ujQ and v^j depending on r. To solve these equations we shall first consider the purely 
algebraic equations (13), (15)-(18). This will lead to the introduction of the susceptibility. 

Introduction of the susceptibility. - By elimination of fy from cqs. (17) and (18) we get 
(iv^-u'^) fQ{r,r',u;\uj)^uj'\^,fx{r,r',u;). (19) 
To obtain fg we have to impose a prescription for the pole at co = co': 



{Q{r,r',u)',uj) 



/2 

V,., 



{lO + i0)2 - UJ'- 



■ fx{r, r', Lo) + s(r, r', to) S{ui - ui') , 



(20) 



with u) + iO a, complex frequency in the upper half-plane and infinitesimally close to the real 
axis. The last term contains an as yet unknown tensor 5{r,r',uj). Substituting eq. (20) in 
(16), and eliminating fp and f/i with the help of cqs. (13) and (15) we get an equality wliicli 
expresses fn in terms of fx and s. If we use (9) to eliminate fn in favor of fs, this equality 
may be rewritten as a linear relation between the coefficients of the polarization density —aX 
and the electric field E: 

-afx{r,r',uj) = eoxir,i^)^E{r,r' ,u) + ^ v^, x(»', w) s(r, r', w) . (21) 

The proportionality constant is the position- and frequency-dependent susceptibility: 

-1 



X{r,uj) = - — 
soP 



1 f 

- '^0 - -2 / 

P Jo 



(W + Z0)2 - OJ'^ 



(22) 



Having solved the algebraic equations of the set (13)-(18) we now turn to the differential 
equation (14). 

General form of the tensorial coefficients. - After elimination of and fp with the help 
of (13) and (15) we get from (14): 



^2 ^2 

Afu{r,r',w) + -^fn{r,r',ui) = --^ [afx{r,r',(j)],j, . 



(23) 



Since fn(r,r',a;) is purely transverse in r, we may write the first term at the left-hand side 
as —V X [V X fn]. Subsequently, we introduce by means of (9), and eliminate fx with the 
help of (21). The result is an inhomogeneous wave equation for f^: 



-V X [V X fE{r,r',(j)] + ^ [1 + x(r,a;)] fE{r,r',(j) = 



pac 



■Vu,x{r,i^)s{r,r',u>). (24) 
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To solve it we introduce the tensorial Green function G of this wave equation, which is defined 
as: ^ 

-V X [V X G(r, r', w)] + ^ [1 + x{r, w)] G(r, r', w) = \5{r- r') . (25) 
In terms of G the solution of (24) reads: 

fB(r,r',a;) = -^ Ur" ^ v'U{r" ,io)Q{ry' ,u;) ■ s{r" y ,u;) . (26) 
J p'a" 

Now that we have found in terms of s it is straightforward to express all other tensorial 

coefficients in s. We first list the results for the coeSicients of the field and polarization 
variables: 

f^(r,r',c.) = -- [fs(r,r',a;)]T , (27) 

fn(r,r',w) = -e^ [^siry , (28) 

fx(r,r',w) = — ^-\v^^x{'r,^)s{r,r',i^)-—x{r,oj)^E{r,r',u), (29) 

fp{r,r',u}) = ^ujv^x{r,u})s{r,r',u)) +ia- [fE{r,r',u))]rj, 

CX CO 

+ 'I^ojx{r,u;)fE{r,r',u;), (30) 
where (26) should be inserted. The coefiicients for the bath variables are 

fY[r,r ,oj ,u) = — 5{r,r ,uj)6{uj -w) - - — — — -2 x(r, w) s(r, r , w) 

7 — 7~?^2 72 X »',w)fj5(r,r ,0; , 31 

pa [uj + iQy - u)' 

fQ(r,r',w',a;) = -ip — fy(r,r', w',a;) . (32) 

The expressions listed here still depend on the tensor s, which is not yet known. It can be 
determined by employing the commutation relations (6). 

Determination of the tensor s{r,r',cu). - The general form (11) should satisfy the com- 
mutation relations (6). The first of these leads upon using eqs. (2) to the following constraint: 

I dr"|fi(r, r", co) ■ fn(r", r', u') - f^r, r" , ^) ■ U{r" , r', u:') 
+f^(r, r", a;) • fp(r", r' , J) - ft,(r, r" , uj) ■ fx{r", r' , J) 

+ dw" ^Ur, r", J', a;) • fQ(r", r', a;", a;') - f^(r, r", J', a;) ■ (r", r', a;", J) | 

= -ih8(ijj -J)\3{r -r'), (33) 

where f^(r,r') equals f*(r',r) with interchanged tensor indices. If the expressions (27)-(32) 
are substituted at the left-hand side, we find contributions with either no factor s, or one or 
two such factors. We first consider the contribution that depends quadratically on s. In the 
r"-integrand in (33) one encounters a frequency integral, which can be rewritten in terms of 
X(r",a;) and x(r", w'). After doing so, one finds that most terms depending quadratically on 
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s drop out. A single term of this type survives. It yields the following contribution to the 
left-hand side of (33): 



2? r ^ 

— S{u-u;') dr" — s^{r,r",Lo)-s{r",r',Lo). (34) 
J P 



The remaining contributions at the left-hand side of (33) can likewise be evaluated by first 
rewriting the frequency integrals in terms of the susceptibility. Subsequently, one finds that 
the integrands contain terms with the transverse part fjj(T", r", lv) of f|;(r, r" , uj). Upon 
adding all such terms one establishes that this transverse part is multiplied by the combination 

[1 + x(r", a;')] fsir", r', u;') + x(r" , J) s(r", r\J) . (35) 

This combination is itself purely transverse in r", as one proves directly from (24). Hence, in 
the integral over r" one may replace the transverse part fjj(r, r", oj) by the full coefficient 

fl;{r,r" ,uj). A similar remark applies to the terms with the transverse part of f£(r", r', w'). 
After these replacements, one arrives at an integral expression, which after use of (24) becomes 
proportional to 



/ 



dr" 



I V" X [V" X ft,(r, r", u;)]| •fs(r", r' , u;')-fUr, r", c.) • | V" x [V" x f^(r", r', u;')]| 



(36) 

A partial integration shows that this integral vanishes. Hence, the only terms that contribute 
to the left-hand side of (33) arc those depending quadratically on s, as given in (34). As a 
consequence, the tensor s has to fulfill the condition: 



/ 



dr" ^ s+ (r, r", uj) ■ s(r", r', w) = ^ 1 5(r - r') . (37) 



Hence, the general form for s is 



s(r,r',u;) = ^^U(r,r',a;). (38) 

The frequency-dependent tensor U(r,r',cj) must be unitary, so that it has the property 
/ dr" U^(r,r",aj) • U(r",r',a;) = \5{r — r'). One may check that the second commutation 
relation in (6) does not lead to new constraints on s. 

Upon inserting (38) and (26) in (27)- (32), and substituting the latter in (11), we obtain 
the diagonalizing operators C{r,Lo) and C^{r,ui). They are determined up to a unitary 
transformation, which leaves both (5) and (6) invariant. The simplest choice for U is a 
local diagonal tensor cxp[i'ip{r, oj)] I 5{r — r'), with an arbitrary phase factor cxp[i'ip{r, lj)]. A 
convenient form for this phase factor is exp[i7/;(r, w)] = ix*('",w)/|x(r,a;)|, as it enables us to 
eliminate v^^ from the tensorial coefficient (26). In fact, from (22) one can prove the identity 

= Q![2pImx(r,u;)/(7reow)]^/2/|x(r,w)|. 

Now that the diagonalizing operators of the model have been found, we can switch to 
the Heisenberg picture and derive explicit expressions for the time-dependent operators rep- 
resenting the electric field and the displacement field. With the above choice for the unitary 
transformation and the associated phase factor we obtain: 



E{r,t) = -i/Jo J dr' dLJOje'^* G{r,r',(j) ■ J{r',uj) + h..c. 



(39) 
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■OC 



D{r,t) 



du! e-''^* w [1 + x{r, w)] G(r, r', u) ■ J{r',iv) 



+i 



r°° 1 

/ dwe""^*- J(r,a;)+h.c., 



(40) 



with a current density J that is proportional to th(^ diagonahzing operator: 



J(r,c.) 



feo Imx(r, cj) 



(41) 



TT 



Similar expressions can be found for the Heisenberg operators that represent other canonical 
variables. 

The expressions (39) and (40) agree with those found before by means of a Laplace- 
transform technique [4]. The electric field E{r,t) at the position r and the time t is given 
by an integral transform, which involves the current density J{r',uj) as a source and the 
Green function G(r,r',a;) propagating a disturbance from r' to r. The displacement field 
D(r,t) is the sum of two terms. First, it contains a contribution that is closely analogous 
to the expression for the electric field E{r,t), with the local permeability 1 + x(r,a;) as an 
extra factor in the integral transform. The second term in the polarization density depends 
on the local current density J{r,u)) only, without an intervening permeability, and without 
propagation effects. For the homogeneous case the expressions (39) and (40) reduce to those 
found in [2]. In a phenomenological quantization scheme [5], forms like (39) and (40) are 
postulated without proof. In that context J{r,u)) is called a noise-current density, which is 
introduced without giving its precise connection to the fundamental dynamical variables of 
the system. In contrast, the present model leads to the specific expression (41) with (11) for 
J{r,uj). 

Conclusion. - By using Fano's diagonalization procedure we have succeeded in finding 
the ladder operators C and C\ which annihilate and create the normal- mode excitations 
of the inhomogeneous damped-polariton model. Knowledge of these fundamental operators 
suffices to determine the full time dependence of the dynamical operators that describe the field 
and the dielectric. The general structure of these time-dependent operators agrees with that 
postulated in a phenomenological quantization scheme. Hence, the present diagonalization 
of the inhomogeneous damped-polariton model gives a justification of the phenomenological 
quantization method for inhomogeneous absorptive dielectrics. 
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